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On the infeasibility of entanglement generation in 
Gaussian quantum systems via classical control 

Hendra I. Nurdin, Ian R. Petersen, and Matthew R. James 



Abstract — This paper uses a system theoretic approach to show that 
classical Unear time invariant controllers cannot generate steady state 
entanglement in a bipartite Gaussian quantum system which is initialized 
in a Gaussian state. The paper also shows that the use of classical 
linear controllers cannot generate entanglement in a finite time from a 
bipartite system initialized in a separable Gaussian state. The approach 
reveals connections between system theoretic concepts and the well known 
physical principle that local operations and classical communications 
cannot generate entangled states starting from separable states. 



I. Introduction 

Entanglement is a unique feature of quantum mechanical systems 
not found in classical systems and is responsible for some of their 
predicted counterintuitive behavior, as exemplified by the famous 
Einstein-Podolsky-Rosen paradox fTl- Entanglement gives rise to 
experimentally verifiable non-classical correlations among measure- 
ment statistics [2| that cannot be explained by the usual classical 
probability models. One well known application of entanglement 
is quantum teleportation, the process of transferring the unknown 
state of one quantum system to another whilst destroying the state 
of the former, without the two quantum systems ever interacting 
directly with one another |3 1. This process is at the heart of quantum 
communication schemes. 

A bipartite quantum system is the composite of two quantum 
systems. The state of such a system will be referred to as a bipartite 
state and is represented by a density operator /qj Suppose that the 
system is composed of a quantum system A with underlying Hilbert 
space Ha and a quantum system B with underlying Hilbert space 
Hb- a state p is said to be separable if it can be decomposed as 
p = 'Y^f.Pkpk ® Pk ^ with p^ and pf being density operators on 
Ha and Hb, respectively, for k = 1,2,.... Here ® denotes the 
tensor product of operators. If a bipartite state is not separable, then 
it is said to be entangled. For a pure state density operator p (i.e., 
tr(p'^) = 1), it can be easily determined if it is separable; e.g., see 
l4l . However, determining the separability of a mixed state bipartite 
density operator p (i.e., tr(p^) < 1), is far from straightforward 
and a complete characterization is only known for certain types 
of bipartite systems, such as for bipartite systems on the finite 
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dimensional Hilbert space C (8) C . In fact, the general problem 
of determining the separability of a given mixed quantum state is 
known to be NP-hard |5|. Another class of bipartite systems for 
which a complete characterization of separability is known is the 
class of bipartite Gaussian systems L6J, UJ, [8J, [9|. These systems 
are commonly encountered in the field of quantum optics. For such 
systems, the underlying Hilbert space is the tensor product of two 
quantum harmonic oscillator Hilbert spaces; e.g., see 1 10 Chapter 
III]. Also, the separability of a state can be completely determined 
from the (symmetrized) covariance matrix of the canonical position 
and momentum operators of the system |6). The class of systems 
considered in this paper is the class of bipartite Gaussian systems. In 
particular, we analyze dynamical bipartite Gaussian quantum systems 
whose covariance matrices evolve in time. Hence, the separability 
or entanglement of these systems also evolves in time. In quantum 
optics, these dynamical bipartite Gaussian systems correspond to a 
class of linear quantum stochastic systems |1I|, |12|, |13| that are 
driven by Gaussian bosonic fields and with a density operator initially 
in a Gaussian state. The dynamics of such systems can be represented 
by linear quantum stochastic differential equations (QSDEs) in the 
canonical position and momentum operators and this makes them 
suitable for a system-theoretic analysis. 

We study the problem of entanglement generation using classical 
finite dimensional (linear time-invariant (LTI) and time varying) con- 
trollers from a system-theoretic point of view. The main contribution 
of the paper is the use of system theoretic arguments and methods 
to show that the application of a classical linear dynamic controller 
cannot induce entanglement in a dynamical bipartite Gaussian system 
which is initially in a separable state. Our result is in agreement 
with the fundamental physical principle that Local Operations and 
Classical Communication (LOCC) cannot generate entanglement 
between initially separable states; e.g., see (H for a proof of this 
result. One motivation for the results of this paper is that they provide 
a starting point for investigating connections between systems theory 
and quantum physical principles. The no-go results for Gaussian 
quantum systems considered here are in a similar spirit to other no 
go results that have previously been obtained in II4I . showing that 
linear modulation of a beam cannot create out-of-loop squeezing, and 
1 15|, showing that neither in-cavity squeezing nor output squeezing 
can be created using linear modulation of the cavity field. 



II. Preliminaries 



A. Notation 



We will use the following notation: i = \/—l, * denotes the adjoint 
of a linear operator as well as the conjugate of a complex number. 
li A= [fljfe] then A* = [a**,], and A'^ = (A*)'^, where ^ denotes 
matrix transposition. ^{A} = {A + A*)/2 and SsjA} = ^,{A - 
yl*). We denote the identity matrix by / whenever its size can be 
inferred from context and use /„ to denote minx n identity matrix. 
Similarly, denotes a matrix with zero entries whose dimensions 
can be determined from context. diag(Mi, M2, . . . , M„) denotes a 
block diagonal matrix with square matrices M\ , M2 , ■ ■ ■ , Af„ on its 
diagonal, and diag„(A/) denotes a block diagonal matrix with the 
square matrix M appearing on its diagonal blocks n times. Also, we 
1 " 



will let J = 



-1 



B. The class of linear quantum stochastic systems 

In this paper, we are concerned with a class of quantum stochastic 
models of open (i.e., quantum systems that can interact with an 
environment) Markov quantum systems that are widely used in the 
area of quantum optics. Such models have been used in the physics 



and mathematical physics literature since at least the 1980's; e.g., see 
fT6l . fTTl . (Tsl, fT^, f2b|. We focus on the special sub-class of linear 
quantum stochastic models (e.g., see 1211 Section 7.2], (20! Section 
6.6], GH Sections 3, 3.4.3, 5.3, Chapters 7 and 10], |22], |23, Section 
5], nil . 1121 . 1131 . 1241 ) that describe the Heisenberg evolution of the 
(canonical) position and momentum operators of several independent 
open quantum harmonic oscillators that are coupled to external 
coherent bosonic fields, such as coherent laser beams. These linear 
stochastic models describe quantum optical devices such as optical 
cavities 1251 Section 5.3.6]t2.1 , Chapter 7], linear quantum amplifiers 
1191 Chapter 7], and finite bandwidth squeezers 1191 Chapter 10]. 
Following nil . 1121 . 1131 , we will refer to this class of models as 
linear quantum stochastic systems. 

Suppose we have n independent quantum harmonic oscillators. 
The jth quantum harmonic oscillator has position and momentum 
operators qj and pj with underlying Hilbert space L^(R); see, e.g., 
1101 Chapter III]. The position and momentum operators satisfy the 
canonical commutation relations [f/j, Pfe] = 2i(5jfe, [gj,gfc] = 0, 
and [pj,Pfe] = 0, where (Jj^ denotes the Kronecker delta and 
[■ , ■] denotes the commutation operator. The quantum harmonic 
oscillators are assumed to be coupled to m external independent 
quantum fields modelled by bosonic annihilation field operators 
Ai{t),A2{t), ■ ■ ■ , Am{t) which are defined on separate Fock spaces 
(over L^(R)) for each field operator 1161 . (l8|. For each annihilation 
field operator Aj{t), there is a corresponding creation field operator 
A'j{t), which is defined on the same Fock space and is the operator 
adjoint of Aj{t). The field operators are adapted quantum stochastic 
processes with forward differentials dAj{t) = Aj{t+dt)—Aj{t) and 
dA'j (t) — Aj {t+dt)—A'j (t) that have the quantum Ito products 1161 , 
fTSl : dAj{t)dAkit)'' = Sjkdt; dA*{t)dAkit) = dAj{t)dAk{t) = 
dA*{t)dAUt) = 0. 

We collect the position and momentum operators in the col- 
umn vector X defined by a; = (qi,pi,q2,p2, ■ ■ ■ ,qn,Pn)^ ■ Note 
that we may write the canonical commutation relations as xx"^ — 
{xx'^)'^ — 2iB with & — diag^(J). We take the composite 
system of n quantum harmonic oscillators to have a quadratic 
Hamiltonian H given hy H = ^x"^ Rx, where 7? is a real 2n x 2n 
symmetric matrix. The quantum harmonic oscillators are coupled 
to the fc-th quantum field via the singular interaction Hamiltonian 
Hk = KLkvUt) - Llr,k{t)) (m, (m, where Lk = KkX 
(with Kk G C^^^") is a linear coupling operator describing the 
linear coupling of the quantum harmonic oscillator position and 
momentum operators to rjk{t). Here rjkit) is a quantum white noise 
process L17J . L19J satisfying the relation Ak{t) ~ /q rik{s)ds. 
We now collect the coupling operators Li, L2, • • ■ , Lm together in 
one linear coupling vector L = (Li, L2, . . . ,-Lm)"^ = Kx, with 
K = [ Ki K2 ■ ■ ■ K^ ]^, and collect the field operators 
together as A{t) = {Ai{t),A2{t),...,Am{t))'^- Then the ;omf 
evolution of the oscillators and the quantum fields is given by a 
unitary adapted process U{t) satisfying the Hudson-Parthasarathy 
QSDE (T61, (13: 

dUit) = ftr{{S - lfdA{t)) + dA{t)'^L - L'^SdA{t) 
-{iH + -L''L)dt\u{t), 



where S £ C™^™ is a complex unitary matrix (i.e., S'' S = SS'' = 
7) called the scattering matrix, and A(i) = [Ajk(t)]j,k=i....,m, with 
Afcj(i) — Ajk{t)*. The processes Ajfc(t) for j,k = l,...,m are 
adapted quantum stochastic processes referred to as gauge processes, 
and the forward differentials dAjk (t) = Ajk {t + dt) — Ajk (t) j, k = 



1, . . . , m have the quantum Ito products: 

dAjkit)dAj,k'it) = Skj'dAjk'it), dAj{t)dAki{t)^5jkdAi{t), 
dAjkdAi{ty = 5kidA*{t), 

with all other remaining cross products between 
Aj(t),Al{t),A,,k'(t) being 0. 

For any adapted processes X{t) and Y{t) satisfying a quantum 
Ito stochastic differential equation, we have the quantum ltd rule 
d{X{t)Y{t)) = X{t)dY{t) + {dX{t))Y{t) + dX{t)dY{t); e.g., 
see 1161 , 1181 . Using the quantum Ito rule and the quantum Ito 
products given above, as well as exploiting the canonical commu- 
tation relations between the operators in x, the Heisenberg evolution 
x{t) = U(tyxU{t) of the canonical operators in the vector x can be 
obtained. Then x(t) satisfies the QSDE (see ||22l, (23l, Cll, (HI): 



dx{t) = d{u{tYxU(t)) 

= Aox{t)dt + Bo 

dY{t) = d{U{tyA{t)U{t)), 
= Cox(t)dt + DodA(t) 



dA(t) 
dA(t)* 



■,x{Q) 



(1) 



with Ao = 29(7? + '=s{K'^K}), Bo = 2ie[ -K'^ S K"^ S* ], 
Co = K, and Do = S. Here, Y{t) = {Yi{t), . . . ,Y,^{t)f = 
U{tyA{t)U{t) is a vector of output fields that results from the 
interaction of the quantum harmonic oscillators and the incoming 
quantum fields A{t). Note that the dynamics of x{t) is linear, and 
Y{t) depends linearly on x{t). We refer to n as the number of degrees 
of freedom of the linear quantum stochastic system. 

In this paper it will be convenient to write the dynamics in 
quadrature form as in |11]: 



dx{t) = Ax{t)dt + Bdw{t); x{Q) = x. 
dy{t) = Cx{t)dt + Ddw{t), 



(2) 



with 

w{t) 

y{t) 



2(K{A(i)},3{A(i)}, 

2(3?{yi(i)},S5{yi(t)},. 



.,?lt{Amm,^{Am{t)}f; 



Here, the real matrices A, B, C, D are in a one to one correspondence 
with the matrices Ao, Bo, Co, Do- Also, the quantity w(f) satisfies 
the Ito relationship dw{t)dw{t)^ = F^dt where 7<V> > 0; see 
1 11|. Furthermore, we define the matrix 5*^, = | [F^^ + 7^*) and 
the differential commutation matrix Tm = | (F,„ — 7*"*). For the 
boson fields that we consider here, T^ is necessarily of the form 
T^ = idiag^(J). The symmetric matrix S^ is then such that F^, > 
and for the Gaussian boson fields that are of interest here this matrix 
reflects the statistics of the field. For instance, Sw ~ I corresponds 
to a vacuum Gaussian boson field. 

Fig. n\ shows an example of a two degree of freedom linear 
quantum stochastic system connected to a classical controller. The 
linear quantum stochastic system consists of two independent optical 
cavities 1251 . 1211 . 1131 denoted by Gi and G2. The two optical 
cavities are connected to the classical controller via a homodyne 
detector (HD) which measures one of the quadratures of the output 
field Yi(t) from Gi, and an electro-optic modulator (MOD) which 
modulates the quantum field .43(f) with the controller output signal 
u{t) and then sends the resulting field .42(f) to G2. 

Remark 1: For the remainder of this paper, we will consider the 
case where n = 2, corresponding to two degree of freedom linear 
quantum stochastic systems, with the quantum harmonic oscillators 
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Fig. 1. The interconnection of two optical cavities Gi and G2 via a classical controller. 



being initialized in a Gaussian stataj Also, for such a linear quan- 
tum stochastic system G, we define a (linear) dynamical bipartite 
Gaussian quantum system (corresponding to G) as the open quantum 
system obtained from G by tracing out (averaging) the bosonic fields. 

III. Simple system-theoretic proof of the Heisenberg 
Uncertainty Principle 

For a system of the form pj, the corresponding symmetric covari- 
ance matrix defined by 

Pit) = ^tr(p(0)(x(i)x(t)^ + {x{t)x{tff)) 

varies with time. Here, p(0) is the initial density operator of the 
overall composite closed system. In this section, we will assume that 
the matrix ^ in J2l is Hurwitz. Then the steady-state symmetrized 
covariance matrix P — limt_s-oo P{t) satisfies the real Lyapunov 
equation (see, e.g., [20' p. 327], 112] Section 4]): 



AP + PA' + BS^B' 



0. 



(3) 



On the other hand, since the commutation relations are preserved, we 
also have that 111.1 : 



Ae + QA' ~ KBT^B' = 0. 



(4) 



Defining the complex Hermitian matrix P = P + iO, we see 
from combining |3b and ([4|( that P satisfies the complex Lyapunov 
equation: AP + PA^ + BF^B"^ = 0, where F„ = S^+T^. Since 
F„ > and A is Hurwitz, it follows that P > 0; e.g., see IZTl . 
Equivalently, in terms of P and O we have that: P + i0 > 0. This 
matrix inequality is a version of the Heisenberg Uncertainty Principle 
that must be satisfied by all Gaussian quantum systems; e.g., see (6l, 
171 for this alternate form of the Heisenberg Uncertainty Principle. 

IV. Classical LTI controllers cannot generate 

steady-state bipartite entanglement in linear Gaussian 

quantum systems 

A. Separability criterion for dynamical bipartite Gaussian systems 

It has been shown in |6| that the separability of a bipartite 
Gaussian density operator p can be completely determined from a 
complex linear matrix inequality (LMI) involving the (symmetrized) 
covariance matrix P — ^tr{p{xx^ + {xx"^)"^)); see also [8J, 191 . 

Lemma 2 (^, ^, ^): A bipartite Gaussian density operator p 
is separable if and only if the corresponding covariance matrix P 
satisfies the LMI P + idiag( J, - J) > 0. 

Note here that without loss of generality, we can assume that x 
has zero mean because the mean of x plays no role in determining 
the separability of the associated density operator. Now, in the case 

^A state with density operator p, is said to be Gaussian if tr(pe'^ ^) = 
gU^m-iA^GA J.JJJ. ^jj y^ g ]g2n^ ^Yisve m e IR2" and G is a real symmetric 
matrix satisfying G + i© > with © as defined previously; see, e.g., 1261 , 
(191,(11. 



of a dynamical bipartite Gaussian quantum system corresponding to 
a linear quantum stochastic system, the covariance matrix can vary 
with time and is given by 

P{t) = ]^tr{p^{t){xx^ + {xx^Y)) 

= \ir{p{Q){x{t)x{tf^{x(i)x{tf)^)), 



(5) 



where p{t) is the density operator at time i > of the overall 
composite closed system, while po (i) is the reduced density operator 
of the two quantum harmonic oscillators at time t obtained by tracing 
out the bosonic fields; e.g., see I.19J . Note that the second equality in 
^ follows from the definition of the partial trace (e.g., see |T8, p. 
102]) since the elements of x are operators on the bipartite quantum 
harmonic oscillator Hilbert space. Also, the final equality in IJSJ 
follows by switching from the Schrodinger picture (in which pit) 
evolves in time) to the Heisenberg picture (in which x{t) evolves in 
time and the overall density operator is fixed as p(0)). Thus, to check 
whether the system is separable at any time f > 0, it is equivalent to 
check if the LMI Pit) + idiag(J, —J) > is satisfied at that time 



B. Separability of dynamical bipartite Gaussian systems coupled via 
a classical LTI controller 

Let Gi = iAi,Bi,Gi,Di) and G2 = (^2,-62,^2,1)2) define 
two linear quantum stochastic systems of the form ^. We form a 
linear quantum stochastic system G = iA, B, G, D) of the form |J2]| 
from Gi and G2, with A = diag(^i,yl2), B = diag(i3i, P2), 
C = diag(Gi,G2), and D — diag(Di,D2)- We then obtain a 
dynamical bipartite Gaussian quantum system corresponding to G 
(see Remark [Tl. The quantum system G is connected to a finite 
dimensional classical controller as shown in Fig.[2]to form a quantum 
feedback control system. In this quantum feedback control system, 
some of the output fields from Gi and G2 are measured and fed to 
the classical controller that processes these measurements linearly to 
produce control signals that are fed back into Gi and/or G2. Here, 
control actuation can be facilitated in two ways: 

1) Modulating the Hamiltonian of Gk by a classical 2x1 signal 
vector ui,k. If the canonical operators of Gk are represented 
by a vector of operators Xk ~ iqk,Pk)'^, this means that the 
quadratic Hamiltonian Hk of Gk is augmented by adding a 
linear (time varying) Hamiltonian term Hi^kit) of the form 
Hi,kit) ~ ui^kit)"^ MkXk, where Mk is a real 2x2 matrix. 
Thus, the total Hamiltonian for Gk becomes Hk + -ff;,fe(i). 
The signal Mi,fe(t) is classical and can depend linearly on the 
classical controller internal variables (i.e., its state) as well as 
the measurement results. This actuation can be implemented in 
different ways, for instance, as described in the Appendix of 
(281. 



2) Modulating (or displacing) an input field of Gk with a classical 
control signal U2,k{t)- This can be implemented by an electro- 
optic modulator. 

In the quantum feedback control system shown in Fig. |2] the vector 
of quantum input fields Ak (t) for the system Gk is partitioned into 
two parts: some of which will be the components of Aki{t) while 
the others will be components of Ak2{t)- Here Aki{t) represents the 
input fields of Gk that will not be modulated by the controller, while 
Ak2{t) represents the input fields that are modulated by the controller. 
Part of the output vector of quantum signals, Yk2{t), of Gk (fc = 1, 2) 
is passed through a network of static optical components (as listed 
in 1131 Section 6.2]) and homodyne detectors (labelled as HDN in 
the diagram) that produces the set of classical measurements signals 
mk{t) which drive the controller. The controller produces two sets 
of classical control signals (k — 1, 2): one set, Uki{t), modulates the 
linear Hamiltonian term Hi^k, and another set, Uk2{t), is modulated 
by a network of (possibly electro-optic) modulators (denoted in the 
diagram by MOD) to produce the quantum signal Ak2{t) as one of 
the input fields into Gk- The signals Vjk(t), j,k = 1,2 are any 
additional quantum noises required for the operation of HDN and 
MOD (they may be suitably absorbed into the definition of wi or 

W2). 



where the real matrices A and B have the special structure: 



YuOl 



m,(i) 



.'^2,(0 



Y:,(l) 




mj(t) 



Fig. 2. Interconnection of Gi and G2 via a classical controller. 

The assumptions that we will use regarding this quantum feedback 
control system are: 

1) The control u{t) = (■Ui,i(t),it2,i(i), wi,2(i), W2,2(i))^ has 
been generated by a finite dimensional linear (time invariant 
or time varying) system. 

2) The quantum signals coming into Gi and G2 come from 
independent sources. Therefore, F^ — dia,g{Fw^,Fw2) and 
hence, 5"^ = diag(5„i, S^a), 7™ = diag(ru,i, T^.^). 

Note that the systems Gi and G2 are not directly connected to one 
another. That is, no output field from Gi is passed directly to G2 
and vice-versa. They are only indirectly connected via the classical 
controller. Note also that the overall closed-loop system is then a 
mixed classical-quantum linear stochastic system as described in |1I| . 

Let z{t) denote the controller internal state which is classical 
in nature and of arbitrary dimension ric- We could also allow the 
classical controller to be driven by an additional classical Wiener 
noise source Wc{t) that is not derived from the measurement signals. 
However, this additional noise may be absorbed into uii or W2; 
see tni for details. Now, let x{t) = {xi{t),X2{t), z{t))'^ and 
w{t) = {wi{t),W2{t))'^ where xi{t) represents the vector of system 
variables for the quantum system Gi and wi (i) represents the vector 
of quantum noise inputs for Gi. Also, X2{t) represents the vector of 
system variables for the quantum system G2 and 102 (i) represents 
the vector of quantum noise inputs for G2. Now, since Gi and G2 
each only interact with the controller, it follows that the dynamics of 
the closed-loop system can be written in the form: 



A = 



All 





Ai3 





A22 


A23 


A31 


A32 


A33 



B = 



Bii 
B22 

B31 ^32 



with j4ii = Ai, A22 ~ A2 , Bii — Bi; B22 = -B2. Our main result 
in this section is the following. 

Theorem 3: Consider any classical LTI controller that is connected 
to the linear quantum system G such that A is Hurwitz in the 
closed loop system ([6]l. Then the resulting closed-loop dynamical 
bipartite Gaussian quantum system is separable at steady state. Thus, 
a classical LTI controller cannot generate an entangled steady state 
from any initial Gaussian state. 

Proof: Since the controller state z(t) is classical, the commuta- 
tion matrix B for x{t) will be degenerate canonical 1 11 Sections II, 

r 01 " 

III, and III C] of the form 9 = 62 

0ncX2 0ncX2 On^X^ 

where Bi = O2 = J. Suppose that the controller state z{t) 
is of arbitrary dimension ric- We have that the closed-loop mixed 
quantum-classical system satisfies the constraint llll Theorem 3.4] 
: Ae + QA^ - xBT^B'^ = 0, with r„, = Amg{T^^,T^^). This 
equation is equivalent to the following: 

Aii0i + 0iAfi - iBiiT^^Bf^ 



A22©2 -l-e2A 

A32e2-iB32r, 



1^227^2^22 



2^22 



®1^31 ~ iSllTiDi^gj^ 

e2Al^ - \B22T^2B32 

—\Bi2'Tw2B32 



(7) 



Multiplying the (2,2), (2,3), (3,2), and (3,3) elements of this 
matrix equation by —1, yields 

Aii0i + 0iAf, - iBiiT^jBf, 



A3i0i-iB3iT„iBfi 



-A2202 - 02Ai2 + 1522^^2 -BJ2 
-A3202+iB32T„2-B?2 



0iA|;-iBiiT„iB|; 

-©2A|2+iS22r™2Sj2 
iB32T„2^32 



= 0. 



(8) 



Letting 6 — diag(6i, — B2, 0„ 



and T„ = diag(T„i, -Tu,2), 



this matrix equality can be written as: AQ + QA — iBT^B — 0. 
We now use the fact that the closed-loop matrix A is Hurwitz. Then, 
as discussed in Section [ni] the symmetrized steady state covariance 
matrix P satisfies: AP + PA^ + BSv^B'^ = 0. Defining P = 
P + iB and F^= S^+f^ 



Lyapunov equation: AP+PA^ + BFu,B'^ 



dx{t) = Ax{t)dt + Bdw{t); x{0) = x, 



(6) 



we have that P satisfies the complex 
0. Recalling that Sw = 
dia.g{Siui, Sw2) and F„ > 0, we note that F^. = S^. + T^^ > 
for i — 1,2. Then, we note that since _F„, = diag{F-uji, F^^), we 
have that F^ > 0; see, e.g., |27|. Therefore, since A is Hurwitz, we 
have that P > 0. Partitioning P according to the partitioning of x{t) 

into its quantum and classical components as P = ^^ ^ 

\_ -r21 "22 

the property P = P + iB > implies that Pn +idiag(Bi, -B2) = 
Pii +idiag( J, — J) > 0. Therefore it follows from Lemmaplthat the 
dynamical bipartite Gaussian quantum system is separable at steady- 



state. Thus, a classical LTI controller cannot generate an entangled 
steady state from any initial Gaussian state. ■ 

Removing the classical controller by defining all its system ma- 
trices to be zero, a special case of Theorem [3] shows that two 
independent and unconnected systems Gi and G2, with Ai and A2 
Hurwitz, which are initially entangled become separable in the steady 
state. 

V. Classical finite dimensional linear controllers 

cannot generate any entanglement in bipartite 

Gaussian quantum systems in finite time 

In the previous section, we have shown that starting from any state, 
separable or entangled, a classical LTI controller cannot generate or 
maintain entanglement at steady state in a linear dynamical bipartite 
Gaussian quantum system. In this section, by a slight modification 
of the arguments of the previous section, we will show that classical 
finite dimensional linear controllers cannot generate bipartite entan- 
glement in a finite time for any initially separable linear dynamical 
bipartite Gaussian quantum system. Moreover, for this finite time 
analysis, we may drop the requirement that the controller is chosen 
so that the closed-loop matrix A is Hurwitz. 

We follow the notation and set up of the last section. Instead 
of considering the steady state covariance matrix P, we now 
consider the symmetrized finite time covariance matrix P{t) = 
|tr(p(0)(a;(t)a;(f)^ + {x{t)x{t)'^)'^)), < t < 00 satisfying the 
Lyapunov differential equation: 

P(t) = AP{t) + P{t)A^ + BS^B'^, P(0) = Po. 

Theorem 4: Suppose that a linear dynamical bipartite Gaussian 
quantum system is initially separable. Then it remains separable for 
all t > under the action of any classical LTI controller. 

Proof: Since the system is initially separable, P(0) +iS > by 
Lemma pi Let P{t) — P(t) + iB. Then following the same lines of 
argument as in the proof of Theorem [3] and by using standard results 
on Lyapunov differential equations, we find that since -P(O) > and 
P™ > (hence also BF^B'^ > 0) that P{t) > 0, for all i > 
regardless of the values of A and B. Similarly partitioning P{t) 
according to the partitioning of x{t) into its quantum and classical 

' Pll(t) Pl2(t) " 
_ Pl2{t) P22{t) _ 

idiag(J, —J) > 0, for all i > 0. This shows that when the bipartite 
Gaussian system is initially in a separable state, then under the action 
of a classical LTI controller it will remain so for all times. ■ 

Remark 5: Note that it is straightforward to extend the proof of 
the above theorem to allow for linear time-varying controllers rather 
than LTI controllers. 

Example 6: Consider the quantum optical system shown in Fig.fTl 
Suppose that the two optical cavities are identical with the partially 
reflecting mirror on each cavity having coupling coefficient 0.01. 
Also, let the position and momentum operators of cavity G, be 
{qi,Pi), and let x = (172,^2,51,^1)^. Let Wii(i) = 2K{A(t)} 
and Wi2{t) = 2Q{Ai{t)}, i = 1,2,3. Then the dynamics of the 
two degree of freedom linear quantum stochastic system (without 
the controller, homodyne detector and modulator attached) is given 
by: 

dx{t) = -0.005x{t)dt ~ 0Ad{w2i{t),W22{t),wii{t),wi2{t))'^ , 

dy°(t) ^ OAx{t)dt + d(w2l{t),W22{t),Wii{t),Wi2it)f . 

Let y°{t) = {yl{t),y°2{t),yl{t),yl{t)f. The amplitude quadrature 
2/3(4) of y°{t) is measured using the homodyne detector and is 
used as the (stochastic) input m\{t) = yo,{t) to a first order LTI 



components as P{t) 



it follows that Pii(t) + 



controUer that produces a two dimensional output signal u{t) 
(ui{t),U2{t))'^ . The dynamics of the controller is: 



dz{t) = Az{t)dt + Bdmi{t), z{0) 



0, 



.(i) = [ a 



cl 



4t), 



where z{t) denotes the state of the controller, and A — —1, B = 
Gi = G2 = 1. The output signal u{t) is passed through an electro- 
optic modulator and sent to the partially reflecting mirror of cavity 
G2. Let x{t) — {x{t), z{t)Y' . We then have that the interconnection 
of the controller with the two cavities via the homodyne detector and 
electro-optic modulator is a mixed quantum-classical system with 
dynamics of the form ^ defined by the matrices 
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and is driven by the noise {wzi(t),w-j,2{t),wi-i_{t),wr2{t)Y' . Sup- 
pose that the bipartite state of the two cavities is in an initially 
entangled bipartite Gaussian state with covariance matrix Pii(O) 
given below in ([9| 



Ai(0) = 



0.5028 





-0.0528 








0.5028 





0.0528 


-0.0528 





0.5028 








0.0528 





0.5028 



(9) 



We take as our measure of entanglement the logarithmic 
negativity En ISD, (D, US). Partitioning Pii(f) into 2 x 

Pll.l(t) Pll,2(t) ' 

Pll,2(t)^ Pll,3(t) 

max(0, - ln( 2v{t))) 



2 blocks as 
by ENJPiiit))' 



-^\JHt) - \/A(t)2-4det(Pn(i)) and A(t) = c 
det(Pii,3(i)) - 2 det(Pii,2(i)). Note that the logarit 



-Biv(Pii(i)) is given 
where u(t) = 



det(Pii,i(t)) + 
garithmic negativity 
is always nonnegative and has a value of zero if and only if the 
state is separable |29J, |4J, otherwise the state is entangled, with a 
higher value of En indicating a higher degree of entanglement. The 
initial value of the logarithmic negativity is _Eiv(Pii(0)) = 0.1054. 
The solid line in Fig. [3] shows that under the action of this classical 
controller, the logarithmic negativity steadily decreases and finally 
goes to zero in a finite time. At this point, the state becomes separable 
and remains so for all future times. If we instead start at an initially 



separable state with covariance matrix Pii(O) as given in 1 10 1 



Ai(0) = 



0.5704 0.0034 0.0562 

0.5704 0.0528 

0.0034 0.6203 0.0499 

0.0562 0.0528 0.0499 0.6203 



(10) 



then the oscillators' joint state remains separable, as shown in the 
dashed line in Fig. [3] 

VI. Conclusions 

By employing system-theoretic arguments and methods, we were 
able to give a systems theory proof of the fact that classical LTI con- 
trollers cannot generate steady state entanglement in linear dynamical 
bipartite Gaussian quantum systems. Furthermore, we also give a 
systems theory proof of the fact that classical linear controllers cannot 
generate entanglement in a dynamical bipartite Gaussian system 




Fig. 3. The evolution of logarithmic negativity EM(Pii{t))ovev time for the system considered in Examplep] The solid line shows the evolution starting 
from an entangled state with covariance matrix Pii(O) in |9j. The dashed line shows the evolution starting from a separable state with covariance matrix 
Pli(O) in (To). 



initially in a separable state. An interesting topic for future research 
is to consider system-tiieoretic analysis of entanglement between the 
continuous-mode output fields. 
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